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Q ■ Abstract 

' We perform the careful study of the gauge copies problem for the direct center 

projection in SU(2) lattice gauge theory. Our results indicate that this gauge is 
not appropriate for the investigation of the center vortices. We also show that the 
point-like objects, nexuses, are important for the confinement dynamics. 

1 Introduction 

The old idea about the role of the center vortices in confinement phenomena |l[ has 
been revived recently with the use of lattice regularization. Both gauge invariant @ and 
gauge dependent J| approaches were developed. The gauge dependent studies were done 
in a particular gauge, named center gauge. Such gauge leaves intact center group local 
gauge invariance. It is believed that gauge dependent P-vortices denned on the lattice 
plaquettes are able to locate thick gauge invariant center vortices and thus provide the 
essential evidence for the center vortex picture of confinement. So far 3 different center 
gauges have been used in practical computations: the indirect center gauge [[| , the direct 
center gauge Q and the Laplacian center gauge ||. It is known that the first two of these 
gauges suffer from gauge copies problem. Many results supporting the above mentioned 
role of P-vortices were obtained in the direct center gauge. Recently the following feature 
of this gauge has been discovered 0: there are gauge copies which correspond to higher 
maxima of the gauge fixing functional F (see below for definition) than usually obtained 
and at the same time these new gauge copies produce P-vortices evidently with no center 
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vortex finding ability since projected Wilson loops have no area law. It has been argued in 
that one can still use direct center gauge to locate center vortices if one uses gauge fixing 
algorithm avoiding "bad" copies of |J. Below we subject this statement to the careful 
check. Another goal of our paper is to investigate properties of recently introduced new 
objects called nexuses || H or center monopoles |10| . One can define nexus in SU(N) 
gauge theory as 3D object formed by iV center vortices meeting at the center, or nexus, 
with the zero (mod N) net flux. We use P-vortices in the center projection to define 
nexuses in SU(2) lattice gauge theory. 



2 Direct center gauge 

Direct center gauge is defined by the maximization of the following functional of the lattice 
gauge field U n ^ Qjj]: 

F (U) = i E (^Trf/^) 2 = ^J2\ ClU-On* + 1) , (1) 

with respect to local gauge transformations, and can be considered as Landau gauge for 
adjoint representation; V is the lattice volume. Condition (|l|) fixes the gauge up to Z(2) 
gauge transformation. Fixed configuration can be decomposed into Z(2) and coset parts: 
Un tfJi = Z ntfl V n ^, where Z n ^ = signTrf/ nM . Plaquettes constructed from Z n ^ field have 
values ±1. Those of them taking values —1 compose the so called P-vortices. P-vortices 
form closed surfaces in 4D space. Some evidence has been collected, that P-vortices in 
the direct center gauge can serve to locate gauge invariant center vortices. It has been 
reported [|j] that projected Wilson loops computed via linking number of the static quarks 
trajectories and P-vortices have area law with the string tension az(2) very close to the 
string tension of the nonabelian theory <Jsu{2)- This fact has been called center dominance. 
Another important observation was that the density of P-vortices scales as a physical 
quantity [|, [TTJ . We inspect these statements using careful gauge fixing procedure. 

The most common method to fix the gauge of the type (|l|) is the relaxation algorithm 
which makes maximization iteratively site by site. The relaxation is made more effective 
with the help of the overrelaxation. It is known that another algorithm - simulated an- 
nealing - is more effective and very useful when gauge copies problem becomes severe |12| . 



Here we do not employ simulated annealing and apply gauge fixing procedure explained 
in details in ref. M. We call it RO (relaxation - overrelaxation) procedure. 

The main problem of the direct center gauge fixing is that the functional F(U) (|l|) has 
many local maxima. We call configurations corresponding to these local maxima gauge 
copies. They are lattice Gribov copies in fact. It is well known that for some gauge con- 
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ditions which are formulated as the maximization of a nonlocal functional (e.g. Landau, 
Coulomb and Maximal Abelain gauges) the gauge dependent quantities depend strongly 
on the local maxima picked up, while to find out the global maximum is impossible. Thus 
it is necessary to approach the global maximum as close as possible. We follow the fol- 



lowing procedure proposed and checked in |12 |: for given configuration we generate N t 
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gauge equivalent copies applying random gauge transformations, and fix the gauge for 
each gauge copy using the RO procedure. After that we compute the gauge dependent 
quantity X on the gauge copy corresponding to the highest maximum of ([[]), F max (N cop ). 
Averaging over statistically independent gauge field configurations and varying N cop we 
obtain the function X(N cop ) and extrapolate it to N cop — > oo limit. This should provide a 
good estimation for X computed on the global maximum unless the algorithm in use does 
not permit to reach the global maximum or its vicinity (the situation we met also in the 
present study). The main difference of the present study from the calculations performed 
earlier is that we use the higher value of the gauge copies (1 < N cop < 20) than it was 
used in refs. [|[ f|, |7|, [Ll| and make careful analysis of N cop dependence. Due to that our 
results differ drastically from those reported previously || f|, |7], [IT| . 

Separately we compute observables using the modified (LRO) gauge fixing procedure 
||: every configuration has been first fixed to Landau gauge, and then the RO algorithm 
for the direct center gauge has been applied. In this case the effect of large number of 
gauge copies, N cop , is not very important, we confirm the results of ref. ||. 

Note that there exists another proposal [|TJj] for the general gauge fixing procedure 
which is free of gauge copies problem. In some particular limit this procedure corresponds 
to the search of the global maximum [[EJ. There is also a class of gauge conditions ||, 
|T3| which do not suffer from the gauge copies problem. 



3 Results 

Our computations have been performed on lattice L 4 = 12 4 for (3 = 2.3, 2.4 and L 4 = 16 4 
for (3 = 2.5. For (3 = 2.3, 2.4 ((3 = 2.5) we study 100 (50) statistically independent gauge 
field configurations. Using the described above gauge fixing procedure we calculate the 
various observables as functions of the number of randomly generated gauge copies N cop 
( 1 < N cop < 20). 

(i) We confirm the conclusion of ref. |J that gauge copies generated via LRO procedure 
have higher maxima of F(U) and thus are closer to the global maximum of F(U). We 
found that F^^ (N cop ) > F^ x (N cop ) for any value of N cop , at any considered value of (3. 

(ii) We find that LRO procedure gives copies with significantly lower density, p, of P- 
vortices than RO procedure. We use the standard definition: p = -^-y Y.n;fi>u(^ ~ Z n ^ v ). 
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Thus gauge copies generated by RO and LRO procedures are indeed different even in the 
limit N cop — > oo. 

(in) The difference between LRO and RO procedure results can be qualitatively ex- 
plained as follows. Fixing the Landau gauge we get the configuration almost without 
P- vortices, the subsequent RO procedure substantially increases the number of P-vortices 
but percolating cluster does not appear. The original gauge field configuration contains 
a lot of P-vortices and the local RO procedure is not able to remove all large (and even 
wrapping) clusters of P-vortices. The field configuration after application of LRO proce- 
dure contains many small P- vortex clusters; the field configuration after application of RO 
procedure contains one large percolating cluster. It seems that this cluster is responsible 
for the area law behavior of the projected Wilson loops (see below). 

(iv) The most important observable is the Z(2)-projected Creutz ratio %(/) which 
we calculate using the procedure suggested in refs. 0, ||. xOO is defined through the 
projected Wilson loops, Wz(2){C) = exp{z7r£(£p, C)}. Here £(£p, C) is the 4D linking 
number of the closed surface, £p, formed by P- vortex and closed loop C. 




Figure 1: The dependence of the Creutz ratios x{I) on the number of gauge copies N ( 



cop 



for (3 = 2.5, L 4 = 16 4 . The error bars are shown for one value of N cop only and they are 
characteristic for the other data points. The dashed line corresponds to the nonabelian 
string tension, asu(2)- 



In Fig.l. we show the dependence of x(I) on N cop for (3 = 2.5. It occurs that this 
dependence is nicely fitted by the function C1+C2/ J N cop . The reason for such dependence 
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Table 1: The comparison of <Jz(2), &su(2) and p for RO gauge fixing center projection. 



N 

1 > cop 


a Z{2)/ 0~SU(2) 


^P/{^SU{2)0 2 ) 




(3 = 2.3 


(3 = 2.4 


(3 = 2.5 


13 = 2.3 


(3 = 2.4 


(3 = 2.5 


3 


0.94(2) 


0.93(2) 


0.98(2) 


1.30(1) 


1.51(1) 


1.74(1) 


20 


0.87(2) 


0.80(2) 


0.83(3) 


1.27(1) 


1.42(1) 


1.61(2) 


oo 


0.82(3) 


0.71(3) 


0.71(3) 


1.24(1) 


1.33(2) 


1.49(2) 



is still to be understood. In Table 1 we give the ratio <yz(2)/ cr su(2) Q- &z(2) is computed 
from x(7) for 3 < / < 4 data at 12 4 lattice and for 3 < I < 6 data at 16 4 lattice. For 
N cop = 3 (number of gauge copies used in ||]) az(2) is close to crsu(2)- But it becomes 
significantly lower for N cop — > oo. Thus RO procedure results strongly depend on N cop . 
It is important that <Jz{2) is 20-30 % lower than asu(2) for N cop — > oo. This implies that 
even if one restricts oneself to RO procedure as it is suggested in |7j , one cannot conclude 
that P-vortices indeed well locate all center vortices. 

(v) For gauge copies generated by LRO procedure we confirm the result of that 
x{I) is zero within statistical errors for any value of N cop . 

(vi) In Table 1 we also show the ratio 2p/a SU ( 2 )a 2 (p is the density of P-vortices). 



As it is claimed in ref. jTT[ in case of the uncorrelated plaquettes carrying P-vortices 2p 
coincides with the dimensionless string tension, <Jsu(2)< l2 - The results presented in Table 1 
show that the density of P-vortices is not proportional to usu(2)0 2 - We have found out 
that for N cop = 3 p is in good agreement with the asymptotic scaling as it was found in 
H . But for N cop — >• oo p deviates from the two loop asymptotic scaling formula. 

(viz) We also investigate the properties of the point like objects, called nexuses. On the 
4D lattice we have the conserved currents of nexuses, defined after the center projection. 
We calculate the phase, Sj, of the Z{2) link variable: Zi = exp(z7rs/), si = 0,1. Then 
we define the plaquette variable op = ds mod 2, (<jp = 0,1). The nexus current (or 
center monopole current [JT0[ ) is then defined as *j = ^5*ap. These currents live on the 
surface of the P- vortex (on the dual 4D lattice) and P- vortex flux goes through positive and 
negative nexuses in alternate order. The important characteristic of the cluster of currents 
is the condensate, C, defined [jOJ as the percolation probability. As it is shown in ref. 
UTOfl the condensate C of the nexus currents is the order parameter for the confinement - 
deconfinement phase transition. We found that C is nonzero for the gauge copies obtained 
via RO procedure (when the projected Wilson loops have the area law). C is zero (in 
the thermodynamic limit L — > oo) for gauge copies obtained using LRO procedure (when 

1 The data for os;y(2) are taken from [|l6) 
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Figure 2: The /3-dependence of the ratio of the nexus condensate, C, to the SU(2) string 
tension in lattice units ITH . 



the projected Wilson loops have no area law). It is interesting that for RO procedure C 
seems to scale as the physical quantity with the dimension (mass) . This is illustrated 
in Fig.2, where we plot the /3-dependence of the ratio C/(a S u(2) a2 ) 2 - Thus these new 
objects might be important degrees of freedom for the description of the nonperturbative 
effects. 

(viii) It is important to perform the same calculations for the indirect center gauge || 
and for the Laplacian center gauge JJ. 

We thank Ph. de Forcrand and T. Kovacs for useful remarks. This study was partially 
supported by grants RFBR 96-15-96740, RFBR 99-01230a, INTAS 96-370 and Monbushu 
grant. 
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